arXiv:hep-ph/0503135vl 14 Mar 2005 


hep-ph/0503135 

MPP-2005-15 


Charged Higgs Boson Pairs at the LHC 

Alexandre Alves^ and Tilman Plehn^* 

^Instituto de Fisica, Universidade de Sdo Paulo, Sdo Paulo, Brasil and 
^Max Planck Institute for Physics, Munich, Germany 

We compute the cross section for pair production of charged Higgs bosons at the LHC and compare the three 
production mechanisms. The bottom-parton scattering process is computed to NLO, and the validity of the 
bottom-parton approach is established in detail. The light-flavor Drell-Yan cross section is evaluated at NLO 
as well. The gluon fusion process through a one-loop amplitude is then compared with these two results. We 
show how a complete sample of events could look, in terms of total cross sections and distributions of the heavy 
final states. 


I. SETTING THE STAGE 

Understanding electroweak symmetry breaking is arguably the biggest challenge in particle physics, and we are confident that 
we will solve it at the LHC. In the Standard Model, a single Higgs doublet gives mass to the up-type and down-type quarks 
and also manifests itself as one scalar Higgs boson. Electroweak precision data indicates that this Higgs boson is light jj]. In 
the supersymmetric extension of the Standard Model we need two Higgs doublets to give mass to the up-type and down-type 
quarks and to cancel anomalies of the fermionic partners of the Higgs bosons. The particle content of a two-Higgs doublet model 
(2HDM) consists of two scalars, a pseudo-scalar and a charged Higgs boson. It has been shown that the LHC is guaranteed to 
hnd one supersymmetric Higgs boson in the weak boson fusion production process with subsequent decay to tau leptons OH. 
However, if we want to test a supersymmetric Higgs sector the charged Higgs boson becomes crucial. While there are many 
ways to add scalars to the Standard Model spectrum and allow them to mix with the Higgs boson, for example radions II, a 
charged Higgs boson with all the appropriate couplings as in the 2HDM is much harder to fake. 


A. Charged Higgs Bosons at the LHC 

Over the years, there have been many studies of charged Higgs bosons at the LHC. The dominant production processes are in 
association with a top quark @>000, in association with a W boson HIS, and charged Higgs pair pro duction |IlI|T^|3 
fldtlLNll^ f^. The most promising decay channels are (depending on the charged Higgs mass) tP iTMll^l^ . bi 
and Wh° l24ll . According to studies by ATLAS m and CMS the process pp —> tH t{Tv) is currently the most 
promising combination, in particular for large values of tan (3. The reason is that the Yukawa coupling of the charged Higgs 
boson to quarks includes a term proportional to mb tan f3, which means that the tH production cross section is enhanced by 
a factor tan^ f3. Unfortunately, the same studies indicate that for mn > rrit the tH production channel fails for tan/3 < 10 
because of the reduced production rate. 

Pair production of charged Higgs bosons is particularly interesting because three different production processes contribute 
to the same hnal state: the usual Drell-Yan production process qq H^H~ through an s-channel Z or a photon 0, the 
loop-induced gluon fusion process gg —> H^H~ IO.H^[T^ . and bottom-parton scattering bb —s- H^H~. The latter two are 
enhanced as cr oc tan"t /3 for large values of tan /3, which can compensate for the loop factor suppression or the small bottom 
parton density. The idea of this paper is twofold. First, in a detailed study (including a next-to-leading order (NLO) calculation 
of the total cross section) of the process bb we establish a reliable prediction for its rate. Second, we show how the 

different production cross sections can be added and how they contribute in the different areas of 2HDM parameter space. In a 
brief addendum we will estimate the supersymmetric QCD corrections to the process bb H^H~ and test the reliability of the 
approximation for the complete set of supersymmetric diagrams by the leading (resummed) tan/3 contributions. 


B. Bottom Parton Picture 

The tan /3 enhancement of the charged Higgs Yukawa coupling, which consists of one term proportional to mb tan (3 and 
another one proportional to mt cot /3, immediately directs us to consider processes with bottom quarks in the initial state. 


Heisenberg Fellow. 
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Figure 1: The dependence of the bottom-exclusive process gg —> bbH^ H~ on the bottom virtuality and transverse momentum. The upper set 
of figures show the variation with charged Fliggs mass, while the lower figures show the distributions for the different contributing subprocesses, 
namely a f-channel top quark and an s-channel Z boson and photon. The s-channel Fliggs contributions are heavily suppressed for this choice 
of parameters. We see that the shoulders of the pT,b distributions extend to values around mir/2 or M/4 if the mass of the final state is 
identified with M = 2m h 13 , 1 ^. 


Conceptually, this leads to the proper description of bottom partons. Lacking a measurement of bottom quarks in the proton, we 
calculate a bottom parton density from the gluon density, with gluons splitting into a bottom quark pair. Because the bottom mass 
is much larger than Aon r), th is splitting can be described perturbatively. It involves a large logarithm, which can be resummed 
using DGLAP evolution 

We illustrate the bottom-parton approach using the process bg —> tH~ 15 1^ . but the arguments given in this section 

are completely general. First, we compute the total cross section for the gluon-induced process (e.g. gg btH~), integrating 
over the final state bottom quark. This class of processes with explicit gluon splitting into bottom quarks, and with a bottom 
jet in the final state, we refer to as bottom-exclusive. In the limit of vanishing mb, the differential cross section with respect 
to the transverse momentum of the bottom quark will diverge like da/dpT,b ~ ^/PT,b- The cross section after integrating over 
the phase space of the outgoing bottom quark will (again for vanishing mb) be approximately proportional to log(p““/p™'"). 
In the presence of a finite bottom mass, mb regularizes the cross section and yields a dependence log(p™“/mb). If the upper 
integration boundary is given by some hard scale in the process, for example the masses of some heavy final state particles 
M = mn + mt ^ mb, the total cross section after integrating over the bottom jet will be proportional to a large logarithm 
log M/mb- This logarithm can be resummed: switching from the bottom-exclusive process gg btH~ to the bottom-parton 
induced (bottom-inclusive) process bg —> tH~ does precisely that. From the above argument it is obvious that the bottom 
parton density as a function of the bottom factorization scale grows roughly like the leading logarithmic contribution log iJiF,b- 

Historically, bottom-parton induced processes overestimated the total cross section if compared to their exclusive counterpart. 
The reason was that the bottom-parton description introduces a new parameter, namely the factorization scale of the bottom 
parton. It corresponds to the maximum transverse momentum of the final state bottom quark which is included in the bottom 
parton density. Since the bottom density is currently perturbatively calculated we can estimate the value of pF,b, because it is 
implicitly part of the definition of the bottom partons: it is the maximum value of pT,b in the bottom-exclusive process (again 
e.g. gg btH~) for which the differential cross section da/dpT,b behaves like ljpT,b and can therefore contribute to the 
bottom parton density — or in terms of the discussion in the previous paragraph: fiF,b = PtV^- Following a naive dimensional 
argument and identifying p,F,b = M implicitly assumes that the asymptotic behavior of the bottom-exclusive process extends 
to = M. If for some reason p^^ < M this will automatically lead to an overestimate of the total bottom-inclusive cross 
section. 

In previous work 1515 1^ we showed that indeed at the LHC it is more appropriate to choose a smaller scale, pF,b = PtT^ ~ 
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M/4, where M is the hard scale in the process, i.e. the mass of the heavy states Mx in the process gg bX or gg bbX. This 
scale choice can be derived in general from the kinematics of the exclusive process at the LHC M . Higher order calculations of 
the processes bb ^ H and bg bH M have shown that this choice of factorization scale is supported by the stability of 
the perturbative series in as- The perturbative behavior of the total cross section bg tH~ also confirms this choice of scales. 
Moreover, a comparison between the NLO results for the exclusive process gg —> bbH El and the one- or two-bottom inclusive 
processes bg —> bH and bb ^ H showed that the different results are now in good agreement Ei- Note that the different cross 
sections should not agree perfectly, as they include contributions at different orders in perturbation theory. Because they differ 
mostly by the fact that the bottom-parton induced processes resum the logarithm log(p““/mf,), the agreement should be better 
for small final state masses, in our case small mn- On the other hand, they should agree within their perturbative error bands, 
because we know the difference between the calculations and we can estimate their respective theoretical errors. 

Going beyond total cross sections, we showed in Ref. H that for the process bg tH~, the distributions of the heavy final 
state particles are not affected by the approximations of the bottom parton picture, namely the assumption that the incoming 
bottom quarks are massless and that they are collinear. These issues will be addressed for bb in sectionHIIl 

Before we analyze bottom-induced charged Higgs pair production, we perform a consistency check. Based on general kine- 
matical arguments, it was shown that for LHC production processes the bottom factorization scale should be chosen around 
M/4 = mrr/2. We repeat the numerical analyses of Refs. and show the results on which this argument is based in 

Fig.lH The two left panels show the intermediate bottom virtuality, the right pair show the transverse momentum of the inter¬ 
mediate or final state bottom. In accordance with earlier results we see that the shoulder limiting the 1/Qb or l/pT,b behavior 
for large momentum transfer sits at smaller values for the transverse momentum than for the virtuality. For different charged 
Higgs masses both shoulders scale with the mass, which is equivalent to the claim that the bottom factorization scale has to be 
proportional to the hard scale g,F,b M = 2mH- The smaller maximum values of Qb and pT,b for larger values of the charged 
Higgs mass is indicative of the limited phase space at the LHC. The value for p^^, up to which we can observe something like 
a ^jpT.b behavior, is clearly smaller than M = 2m h and in agreement with the general estimate pp^b ~ M/4 = mrr/2. Also 
in Fig. n we see that the shoulder for the s-channel process extends to larger values of Qb an d vt b by about a factor of 1.5. 
This can be understood from the approximate determination of the factorization scale in Ref. 11^ : the position of the shoul¬ 
ders in the Qb and the pp.b curves can be approximately computed relying only on the parton x behavior of the gluon density. 
However, the results ~ M/2 and ~ M/4 assume that the heavy final state is produced at threshold. The difference 

between the s-channel and f-channel processes is that f-channel Higgs pair production can proceed through an S wave, which 
cross section increases proportional to the Higgs velocity (3 El at threshold, while the s-channel process increases like (3^ and 
therefore shows a delayed onset. The s-channel process is strongly suppressed at threshold and produces the bulk of the Higgs 
pairs at a larger partonic center-of-mass energy (we show the corresponding pt.h distributions in sections miandllVl. Thus, 
for the s-channel process, the typical final state invariant mass is larger, which means effectively M > 2m h- This is what we 
observe for the gauge boson exchange process in Fig.^ We also notice that the determination of the bottom factorization scale 
has a theoretical error, because we have to approximate the physical curves with a box-shaped object. Effects like threshold 
suppression of the s-channel diagrams are covered by the theoretical error, and we discuss the different scale variations in detail 
in sectionim 


Organization of the Paper and Conventions 

The paper is organized as follows: in this section we reviewed the status of charged Higgs searches at the LHC and the 
status of the bottom-parton approach. In section Gil we compute the bottom-parton induced cross section to NLO. We show 
how the description is perturbatively stable with respect to QCD corrections. In sectionGIIlwe use the NLO predictions for the 
distributions of the heavy final states to prove the validity of the bottom-parton approach. The results from these two sections 
allow us to quantitatively compare the different production processes for charged Higgs pairs at the LHC in section GYI This 
comparison covers total rates as well as differential cross sections. As an addendum we compute the full SUSY-QCD corrections 
to the bottom-parton scattering process in sectionIVl 

Before we proceed with the discussion of the results we briefly give all conventions used for our calculations. All bottom- 
parton induced cross sections are computed with zero bottom mass. When we check our results using final state bottom quarks 
we use the on-shell mass mb = 4.6 GeV, or for testing purposes a second value 0.46 GeV. The top quark mass is always chosen 
as mt = 178 GeV. In contrast, the bottom and top Yukawa couplings are defined in the MS scheme, i.e. numerically we use 
mb{pR) and mt{pR). Unless explicitly mentioned we use the renormalization scale pR = mu and the bottom factorization 
scale pF^b = fnH The parton densities are CTEQ6 El, consistently chosen to be the LO or NLO sets LI or M; as is also 
calculated at equivalent order. 
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Figure 2: LO and NLO total cross section for the process bb H~ at the LFIC. The K factor is defined consistently as crNLo/crLO- The 

parameters in the 2HDM Higgs sector are chosen consistently for each value of tan j3 and charged Higgs mass niH- In the lower two rows the 
contributions from s-channel Higgs and 7 , Z exchange and f-channel top quark exchange are shown. In those figures the NLO curve is always 
the upper one in each pair. 


II. BOTTOM-INDUCED PRODUCTION PROCESS 

Before we can discuss the different charged Higgs pair production processes at the LHC (in section HVI l. we have to find a way 
to reliably predict the bottom-induced production rate bb H^H~. This shortcoming was noticed in an earlier paper and 

with the recent progress in understanding bottom-parton induced processes we can now fill this gap. The key ingredient for the 
correct evaluation of production processes of the kind bb —> H^H~ is the choice of bottom parton factorization scale, which 
we reviewed in detail in section irm Next, we will compute the NLO corrections (for now in the 2HDM, i.e. without gluino 
contributions, which are added in sectionlVl. 

We calculate the 0{as) corrections to the process pp/bb H'^H~ + X following two different approaches. First, we use 
a one-cutoff technique to calculate the total cross sections. These results will be made publicly available as part of an extension 
to the Prospino2 package 1^ . This approach has the advantage of being numerically fast and reliable. In a second step we 
repeat the calculation with two cutoffs for the soft (Sg — 10“"^) and collinear (6c = ^s/100) divergences 1^ . This allows us 
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Figure 3: The rapidity and transverse momentum distributions at LO and NLO. The results are shown for two different values of tan/3, so 
illustrate the behavior of the s-channel and the t-channel diagrams. We include a set of NLO curves for the factorization scale = 2fi% = 
niH to illustrate the remaining theoretical uncertainty. 


to calculate distributions for the heavy final state particles. These distributions are a crucial check of the validity of the bottom- 
parton approximation, which we use to improve the prediction of the total cross section. Moreover, we will be able to compare 
the differential cross sections for the different charged Higgs pair production mechanisms in section HVI 

NLO corrections: We present the results of our NLO calculation in Fig. |2] The total cross sections are plotted as a function 
of the charged Higgs mass for four different values of tan/3. The 2HDM parameters are chosen consistently 1^ . Before 
discussing the higher-order corrections we note that indeed the leading order (LO) cross section increases like tan"' /3 for large 
values of tan (5. The K factor in the second panel is defined as unlo/<7lo- The size of the NLO correction is strongly dependent 
on the Higgs mass and on the value of tan /3. To understand this we separate the three different contributions to the production 
rate. 

The bottom-induced Drell-Yan contribution proceeds through an s-channel photon or Z boson. The coupling to the pho¬ 
ton is independent of tan /3 and the dependence of the Z coupling is small for large enough values of the pseudoscalar and 
charged Higgs masses (the decoupling limit). The NLO corrections for the bottom-induced Drell-Yan process vary around 
K = 1.60 • • ■ 1.55 for Higgs masses from 160 GeV to 500 GeV. This correction is somewhat larger than for the light-flavor 
Drell-Yan process, which is due to the x dependence of the bottom parton densities and can also be observed in the perturbative 
behavior of the Standard Model bb ^ H process iH. 

A second set of s-channel production processes is made possible by the incoming bottom quarks and their finite Yukawa 
coupling; we can couple the charged Higgs pair to s-channel neutral CP-even Higgs bosons CP-odd pseudoscalar 

exchange is forbidden by the CP symmetry of the final state. The Higgs self-couplings hPH^H~ and are strongly 

dependent on the Higgs sector d' parameters, but the relative size of the Higgs contribution compared to the Drell-Yan process 
is basically negligible; the suppression varies from 10“^ for small Higgs masses below 200 GeV to lO”"^ for Higgs masses 
around 500 GeV. While it is possible that there are regions of 2HDM parameter space where this hierarchy is less pronounced, 
the impact of the s-channel Higgs exchange is small for our analysis. We note, however, that the NLO corrections to the s- 
channel Higgs exchange are smaller than for the Drell-Yan process, namely Ar=1.15--1.20 (for charged Higgs masses from 
160 GeV to 500 GeV). This is in part an effect of the MS renormalization of the bottom Yukawa coupling. The combination of 
the two s-channel processes receives a fairly constant NLO correction, from K = 1.58 for small Higgs masses and tan /3 = 10 
to K = 1.55 for large values of tan (3 or large Higgs masses. 

For the total rate the s-channel production processes compete with f-channel top quark exchange. The btH'^ coupling consists 
of two contributions, yf,tan/3 and yt cot /3. Over the entire parameter space tan/3 > 10, the contribution dominates over 
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the contribution by at least one order of magnitude. We note that because of the different running of the bottom and top 
Yukawa couplings the K factor of the yf contribution is considerably smaller than for the yf contribution @|. This leads to a 
tan /3 dependence of the NLO corrections to the t-channel process. For small Higgs masses the size of the correction varies as 
K = 1.06 ■ • ■ 1.36 (for tan/3 = 10 to 50). For larger values of the Higgs masses (and thereby the renormalization scale) the 
behavior of K becomes flat. 

Combining the s-channel and f-channel contributions, we see that they interfere destructively. For small Higgs masses the LO 
rates are almost identical in size for tan /3 ^ 20. Because the two classes of diagrams have very different NLO cotTections this 
leads to huge relative QCD corrections to the total rate. We observe this in Fig. |2l where the K factors are fairly constant for 
the s-channel dominated regime tan (j < 10, and for the t-channel dominated regime tan /3 > 30. The curve for tan /3 = 30 
shows a small remaining increase only for very small Higgs masses. The slight difference between the large uih behavior of the 
tan /3 = 30 and tan /3 = 50 curves is also due to the interference between the s- and t-channel subprocesses. 

We emphasize that in this paper we show results only for large charged Higgs masses mn > rrit. For small Higgs masses 
there is a technical complication, because the production process bg tH~ and its charge conjugate will be part of the NLO 
corrections. We subtract the cotTesponding divergence consistently in the small width approximation, which is also used in 
Refs, The public computer code in the extension to Prospino2 will, however, include this region of parameter space. 

Phenomenologically it is less interesting, because charged Higgs bosons will be produced in large numbers in anomalous top 
quark decays. Matching with the off-shell process bg tH~ is studied in detail in Ref. 0. 

Distributions: The effect of NLO corrections on the distributions of the heavy hnal state Higgs bosons is shown in Fig. |3 
The case tan/3 = 50 corresponds to the parameter point SPS4 1^ . and the parameter choice tan/3 = 10 is the same in 
all parameters except for tan /3. According to Fig. |2] t-channel top quark exchange is effectively switched off this way. All 
that is left is s-channel gauge boson exchange. For tan (3 = 50, the rate is completely dominated by t-channel top exchange: 
the relative size of the two contributions, (Ts_ only/o’*, only, varies from 2.0% to 1.0% at LO and from 2.5% to 1.5% at NLO 
(with a charged Higgs mass between 200 GeV and 500 GeV). While the shift in this fraction is sizable, as we would expect 
from the different K factors described above, any distribution for tan /3 = 50 will be cornpletely independent of the s-channel 
contributions. Similar to higher order corrections to supersymmetric particle production we indeed see that a constant K 
factor in the {pr.H — Uh) plane is a reasonable approximation within any error estimate. 

The situation is slightly different for tan /3 = 10. Here the peak in the pt.h distribution shifts by 20 GeV when we switch 
to the slightly harder NLO description. The interference between the s channel and the t channel subprocesses plays no role — 
O'*, only/0’s, only decreases from typically 0.25 at LO to NLO values around 0.18, which we checked is not sufficient to impact the 
distributions visibly. Instead, the shift towards a slightly harder spectrum is a result of the three-particle phase space at NLO. The 
extra jet balances the slightly harder charged Higgs bosons in the final state. The same effect actually appears for tan j3 = 50, 
but it is less pronounced. To estimate how the shift in the pt,h distribution compares with the theoretical uncertainty, and to 
understand the difference in the size of the shift between s-channel and f-channel diagrams, we change the bottom factorization 
scale to twice its usual value: pp = mn- From section ITSl we know that this will make the s-channel (P-wave) subprocess 
more similar to f-channel (S-wave) top quark exchange. Indeed, in Fig.|3we see that the LO curve shows no visible effect, and 
that the NLO and the LO distributions are now almost identical. 

As we pointed out in section irm the determination of the bottom factorization scale shows an inherent theoretical error, 
because it is not possible to tell where exactly the asymptotic behavior of the bottom-exclusive processes is no longer valid. The 
observed scale-dependent shift between the LO and NLO distributions reflects this remaining theoretical uncertainty, and shows 
that it is perturbatively under control. 

Scale dependence: Even though we have argued before that the bottom factorization scale should be chosen around mrr/2, 
the renormalization and factorization scales are nevertheless free parameters introduced by the perturbative expansion of the 
hadronic production cross sections. The bottom factorization scale plays a special role, because it is not actually extracted from 
data. Instead, it is calculated from the measured gluon density in a perturbative regime. However, while we can claim what 
the (roughly) appropriate values for ppfi should be, it is also obvious from Fig. [TJthat a serious theoretical etTor margin should 
be assigned to the determination of the appropriate scale. We therefore proceed to check the factorization and renormalization 
scale dependence of the cross section pp —> H'^H~ at LO and NLO. Because the scales are artifacts of perturbation theory, we 
expect the scale dependence to be a measure of the theoretical error on the total cross section, and we expect it to flatten after 
we include the NLO contributions. 

The scale dependences of the LO and NLO cross sections are shown in Fig. 0 For the two dimensional plots both scales 
are varied up to factors 1/10 and 10, respectively, around their central values p^p = mpil^ and = m^f. The central value 
of the factorization scale was motivated in section HE and the central value of the renormalization scale is motivated by the 
production of heavy supersymmetric particles at the LHC M and by the production of Higgs bosons from bottom-parton 
initial states IHIa I^ . Note that the renormalization scale covers the running top and bottom Yukawa couplings at LO, whereas 
the strong coupling as only enters at NLO. 
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Figure 4: Factorization and renormalization scale dependence of the total cross section. The scales are varied around the central values 
f/p = niH /2 and = niH- In the second row the contributions to the NLO rate are separated by their different initial state. In the third row 
the scale, which is not varied, is fixed to its respective value //’. The LO curves are the dashed ones and the NLO curves are the solid lines. 
For the charged Higgs mass we use our central parameter point SPS4 with tan /3 = 50 and tuh = 342 GeV. 


From the two-dimensional behavior shown in Fig. |3 we see that the NLO rate has a considerably flatter scale dependence 
than its LO counterpart. There is, however, one exception to this rule: for very small iip and very small fin the perturbative 
expansion of the process bb —> H'^H~ becomes unstable. This is the regime where the bottom-parton description breaks down: 
a small bottom factorization scale means that hardly any phase space of the gluon splitting in the bottom-exclusive process is 
resummed and very little of the bottom-exclusive production rate is included in the bottom parton density. This means that the 
fcg-induced process which contributes to the NLO rate will dominate over the 66 -induced LO process. This effect becomes more 
pronounced if the perturbative expansion parameter which suppressed the 65 -induced channel as compared the bh-induced Born 
process is large, i.e. for a small renormalization scale which increases the value of asifin). The same feature can be seen in the 
associated production of a charged Higgs with a top quark 01 ■ 

If we stay away from this extreme (and unphysical) scale choice, we can attempt to estimate the theoretical error of the LO 
and NLO total cross sections for pp/bb —> H'^H~ + X. A reasonable range of scales could be a band from /i°/3 to 3/i° around 
the central values. First, we see that the LO and NLO error bands for both scales would overlap in Fig|3 This is not crucial, 
but it is certainly a welcome feature, which indicates that there are no large constant terms entering the cross section calculation 
at the NLO level. Moreover, we see that the factorization scale dependence covers a little less than a 20% error band, and the 
renormalization scale points to an error slightly above 10%. Even though the two scale dependences cancel each other if we 
identify fin/pPn ~ Ff/Mf prefer to conservatively estimate the overall error by adding the two errors in square. This leads 
us to an estimate of the theoretical error (due to higher order corrections) for the NLO total cross section prediction of around 
25%. We emphasize that this estimate cannot be a precise error band with a statistical meaning, as for example statistical or 
systematic experimental errors are. The range of scales we choose to cover is not fixed by first principles, and it is by no means 
guaranteed that no large finite terms come in at the NNLO level. However, we know that we have resummed the leading large 
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Figure 5: The rapidity and transverse momentum distributions of the final state charged Higgs boson for different processes. The two-bottom 
and one-bottom exclusive results (with the initial states gg and bg) are given for the physical mass mt, = 4.6 GeV and for a hypothetical value 
of 0.46 GeV. The NLO results for the 66-initiated process is shown for massless bottom quarks. 


logarithms into the bottom parton densities, and we see a well behaved perturbative behavior going from the LO to the NLO 
cross section calculation. 


III. BOTTOM PARTON APPROXIMATION 

The bottom-parton approach which we use to compute the double-bottom inclusive process bb H^H~ resums large 
logarithms and improves the prediction of the bottom-inclusive cross sections, but it also implicitly makes two approximations. 
One is that incoming bottoms are massless, the other one is that the bottom jets over which we integrate are collinear. For total 
cross sections these approximations are well founded; bottom mass effects are suppressed by powers of rrib/mH, so they are 
well below any errors we quote on the NLO total cross section calculation. The kinematical configuration of the final state does 
not enter the computation of the total cross section, which means that we do not expect any non-negligible finite-pT,b effects for 
the total cross section calculation using bottom partons. 

In section irm we showed that the bottom-parton approach works well for total cross sections. In some sense this is by 
construction, because we estimate the value of the only new parameter, the bottom factorization scale, such that the total cross 
section prediction is consistent with the definition of the bottom parton densities as they are computed from the bottom-exclusive 
process. The only serious test of the bottom-parton approach is the check that our estimate of the factorization scale indeed 
predicts oc M, where M is the hard scale in the process, usually chosen as the sum of final state masses. Differential cross 
sections, in contrast, are a non-trivial check of our approach. Phenomenologically it is crucial that the distributions of the heavy 
final states are not significantly affected by the two approximations mentioned above (here, significantly means with a larger 
effect than the perturbative error of the NLO calculation, which we estimated to be around 25% for the total cross section). 

Massless bottom quark approximation; The first approximation we have to check is the assumption that the bottom partons 
and bottom jets involved are assumed to be massless. This is easy to estimate, because we can compute the one-bottom or 
two-bottom exclusive processes at tree level with different bottom masses. There is a technical complication in the one-bottom 
exclusive processes bg —> bH^H~, where we have to choose the incoming bottom as massless and the outgoing bottom as 
massive. Because of SU(3) gauge invariance we cannot couple a massless bottom and a massive bottom to a gluon. However, in 
our s and f-channel diagrams we are not sensitive to U(l) or SU(2) gauge invariance, so we can switch from massless incoming 
bottoms to massive outgoing bottoms at the bhZ and 667 vertices. 
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Figure 6: The rapidity and transverse momentum distributions for the final state charged Fliggs boson for different choices of the factorization 
scale. The small scale choice = 20 GeV interpolates into the region where the formally-NLO subprocess bg bH^H~ dominates the 
NLO rate. 


We compute all distributions for the physical on-shell mass mb = 4.6 GeV and for a hypothetical 0.46 GeV. Apart from the 
NLO predictions for the bottom-inclusive process bb H^H~, the normalized distributions for the rapidity and the transverse 
momentum of the final state Higgs bosons are computed for the processes gg —> bbH^H~^ and bg bH^H~. These five 
curves are shown in Fig.|3 Because the s-channel neutral Higgs bosons do not contribute visibly, one can think of the different 
values of tan /3 as a rescaling of the independent Yukawa couplings; the parameter choice tan /3 = 10 switches off f-channel top 
quark exchange and leaves us with s-channel gauge boson exchange. 

In Fig. 1^ we see that the bands defined by varying the bottom mass for a given process are considerably smaller than the 
difference between the different initial states, i.e. the difference between the bottom-induced and gluon-induced processes. We 
also see that for smaller values of mb the Higgs distributions become slightly softer. While naively the finite phase space would 
move the pt,h distributions to a softer regime when we increase the bottom mass, the opposite happens: a smaller bottom mass 
in the final state jet cuts off the pT^t distributions at smaller values of pT,b- Because for both bottom-exclusive processes the 
charged Higgs bosons balance the final state bottom jet in the transverse plane, the pr spectrum for all particles becomes softer 
when the bottom jets become softer, which is the case for a smaller cutoff value p™'" ~ mt- However, this effect is already part 
of the collinear bottom quark approximation, which we consider next. Therefore, at this point we conclude that bottom quark 
mass effects are negligible compared to the collinear bottom quark approximation. 

Collinear approximation: From the above considerations we know that the dominant error is the assumption that the final state 
bottom jets are collinear and are resummed into the bottom parton densities. We can test the validity of this assumption using 
the NLO differential cross sections for the process bb part of the NLO corrections are real jet emission diagrams, 

either real gluon emission from the Born process or the bottom quark emission bg —> bH'^H~, which we usually refer to as 
crossed diagrams. In Fig. 0 we saw (and discussed in detail) that, in the limit of pp ^ mb, the perturbative expansion around 
the LO process bb H'^H~ fails. The bottom parton density becomes small, because it is defined to include the splitting of 
the gluon into a bottom pair up to pp.b < Pf- With the lower cutoff this means that less gluon splitting is included in 

the bottom parton density and more of it has to be considered explicitly, which in turn means accounting for the gluon-initiated 
process bg —> bH'^H~, and eventually gg bbH^H~. As we have seen, the K factor for the process hh H^H~ blows up 
in this limit, because already the NLO contribution bg —> bH'^H~ becomes dominant. 

This behavior illustrates how the variation of pp interpolates consistently between two regimes: for large pp the LO process 
bb H^H~ is dominant. In fact, for pp b = m/f/2 it is designed precisely to minimize the effect of the NLO contribution 
bg —> bH'^H~ 1^ 1^ . For smaller pp the formerly NLO contribution bg bH^H~ will become dominant. The NLO 
calculation of hb —> H^H~ combines these two processes consistently (using the zero bottom quark mass approximation. 
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Figure 7: Production cross sections for pair production of charged Higgs bosons. The curves for the bottom-induced process (sum of s- and 
f-channel) is labeled bb. The Drell-Yan process with only light-flavor incoming quarks is labeled DY, and the gluon fusion process through 
the one-loop amplitude is labeled gg. The cross sections are given for different values of tan /3. The bottom-induced production rate and the 
Drell-Yan production rates are computed at NLO. 


which we tested above) @ 11 . We should note that a NNLO calculation including the gg hhH^H subprocess would of 
course be perfectly suited for this test, but goes beyond the scope of our paper. Moreover, we emphasize that if the bottom-parton 
approach using the collinear bottom quark approximation should not work well, this behavior will already show up in the NLO 
interpolation. 

In Fig. 1^ we show the numerical results of this comparison. As usual, we test the assumption for two values of tan/3 or in 
other words independently for the s-channel and (-channel subprocesses. All shifts from LO to NLO curves were discussed 
in section |n] To test the bottom-parton approximation we now change the factorization scale to fip = 20 GeV. We see that 
for both the s-channel and (-channel processes the pr distributions become slightly harder when we compare the leading order 
and the low-scale curves. The now-dominant three particle kinematics is responsible for this shift. In an earlier work we 
showed that part of this effect is also due to the hardening gluon spectrum when we move to smaller values of fip = fi'F,g 
in the bottom-gluon initial state 0]. We note, however, that the shifts in the pt,h and yn distributions are well within the 
perturbative uncertainty. As we discussed in section |II| the NLO distributions for the central factorization scale pp = ruH 
have to interpolate between the LO and the low-scale curves, because they are affected by the two- and three-particle kinematics. 
It is a curious coincidence that the NLO rapidity distributions for this scale choice are much more similar between the s- and 
( channels than any of the other cases. We therefore conclude that the collinear approximation in the bottom-parton picture is 
indeed valid for the distributions of the (heavy) final states. We have made similar checks for larger Higgs masses, and both 
the zero-mass and collinear approximations improve in accuracy. This is in agreement of what one would expect from the 
bottom-parton approach S. 


IV. COMBINATION OF PRODUCTION PROCESSES 

Charged Higgs boson pairs can be produced in three different ways at the LHC |In|. Theoretically the most straight-forward 
production process is light-flavor quark scattering via an s-channel photon and Z lllll . As described in section |II| the charged 
Higgs couplings are nearly independent of tan /3, as long as the pseudoscalar Higgs mass and the charged Higgs mass are 
large enough. We calculate the NLO production cross section in complete analogy to the production of two supersymmetric 
sleptons, a production process which is part of the publicly available computer program Prospino2 m. The factorization and 
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Figure 8: The rapidity and transverse momentum distributions of the final state charged Higgs boson for the three different production mecha¬ 
nisms. The curves for the bottom-induced process are labeled bb. The Drell-Yan process with only light-flavor incoming quarks is labeled qq, 
and the gluon fusion process through the one-loop amplitude is labeled gg. The bottom-induced production rate and the Drell-Yan production 
rates are computed at NLO. 


renormalization scales are chosen as niH- In Fig.0we see that for charged Higgs masses between 160 and 500 GeV the NLO 
rate Drell-Yan process varies between 58 fb and 0.23 fb. This rate is computed without any cuts on the final state particles. 
The size of the NLO corrections for this process varies from K = 1.27 to K = 1.17 (for Higgs masses from 160 GeV to 
500 GeV). The remaining theoretical error on the NLO cross section can be estimated to be less than 25% l30ll . For small values 
of tan P <20 this production process is the dominant channel at the LHC. 

Similar to the light-flavor Drell-Yan process, charged Higgs pairs can also be produced in bottom-parton scattering. We 
analyzed this production process in detail in this paper and can now quote a reliable estimate of the production cross section. We 
And that the s-channel and the f-channel subprocesses interfere destructively. The LO total cross section is enhanced by tan"^ /3 
for large values of tan/3. In Fig. 0 we see that even for tan/3 = 50 the NLO cross section is still below 10 fb. In this region 
of parameter space the f-channel subprocess dominates completely. Compared to earlier analyses we And that we have to use 
running MS Yukawa couplings and a low bottom factorization scale, which suppresses the cross section to a level where it will 
always stay below the light-flavor induced rate as well as below the gluon-induced process discussed in the next paragraph. The 
remaining theoretical error on this total cross section we estimated in section|n|to be on the order of 25%. 

For large values of tan/3, the one-loop process gg H^H~ cross section also increases like tan^/3. In the 2HDM this 
process proceeds through a mixed top-bottom loop or though an s-channel neutral Higgs boson. Even though our estimate 
of the gluon fusion rate is less precise than for the other two channels, we use our LO estimate to compare the power of the 
different channels. The theoretical uncertainty we would attach to this LO cross section prediction is probably around a factor of 
two. However, we believe that quoting the LO rate is conservative: we expect the NLO corrections to be dominated by initial- 
state radiation and therefore similar to the neutral Higgs pair production counterpart As is seen in numerous processes, 
this universal initial state radiation leads to large positive corrections. We also note that this process is not evaluated with MS 
Yukawa couplings, but with their on-shell definition. For loop-induced processes of this type, this renormalization scheme leads 
to a well-behaved perturbative series 1?^ . The gluon fusion rate has a minimum around tan/3 = 7, and in Fig. 0 we see that 
it passes the (constant) Drell-Yan process at the LHC for tan /3 > 40. As stated before, the bottom-parton scattering rate is 
always smaller than the loop-induced gluon fusion rate for production of charged Higgs boson pairs at the LHC. 

Naively adding the rates for the three production mechanisms for charged Higgs boson pair production is consistent, and 
it tells us which process contributes what fraction of the events to the total cross section. However, it will not be the entire 
story. Detector analyses require acceptance cuts on the Higgs decay products, and these cuts will also probe the distributions of 
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Table I: Supersymmetric corrections to the production cross section bb from the resummed At, corrections and explicit remaining 

supersymmetric loop diagrams. The supersymmetric parameter points are chosen according to the benchmarks in Ref. 01. All masses are 
given in units of GeV. The percentage changes are defined with respect to 2HDM rates at NLO. 


the charged Higgs bosons. Therefore, we have to compare the transverse momentum and rapidity distributions of the charged 
Higgs bosons in the three different production mechanisms. As before, we use the NLO results for the bottom-parton scattering 
process and the Drell-Yan process. The comparison is shown in Fig.|3 As usual, we show the results for two different values of 
tan (3, which distinguish the s-channel (tan /3 = 10) and the t-channel (tan /3 = 50) contributions. In contrast to bottom-parton 
scattering, we see that the distributions for the Drell-Yan and gluon fusion processes do not depend on tan /3. For Drell-Yan 
s-channel photon and Z exchange this reflects the fact that the entire scattering matrix element is to a good approximation 
independent of tan/3 (for large enough mn)- In gluon fusion, the steep x behavior of the gluon parton distributions pushes 
the produced heavy particles close to threshold, which means the two charged Higgs bosons appear with very little transverse 
momentum. Again, this feature is independent of tan/3. The bottom-parton scattering spectrum varies strongly with the relative 
weight of the s-channel and the t-channel diagrams. At threshold, f-channel top quark exchange (^-wave) increases rapidly 
with /3, much faster than s-channel gauge boson exchange (P-wave), which shows a delayed increase, as /3^. The transverse 
momentum spectrum becomes softer the more the f-channel contributes, i.e. for larger values of tan /3. This effect is particularly 
strong because the bottom parton densities are derived from the gluon density and therefore are steep in x, pushing the partonic 
center-of-mass energy to the smallest possible values. 


V. SUPERSYMMETRIC CONTRIBUTION TO BOTTOM-INDUCED PROCESS 

Diagrams with a gluino instead of a gluon contribute to the NLO rate for the process pp H^H~ . However, for massive 
gluinos only virtual diagrams appear, in which the gluon and the quark lines are replaced by gluinos and squarks. Usually, the 
SUSY contributions are suppressed by the large supersymmetric partner masses. However, there is a well-known exception to 
this rule, where the power suppression by large sbottom masses can be compensated by the gluino mass and the Higgsino mass 
parameter in the numerator. Sbottom-gluino corrections to the external (massless) bottom legs can involve a left-right mixing 
for the sbottom propagator. Because the off-diagonal matrix element in the sbottom mass matrix has the form rub {Ab — p tan /3) 
the left-right mixing is enhanced by a factor tan /3. These potentially large corrections will automatically become important for 
charged Higgs searches whenever we rely on a large production rate with values tan /3 > 30 in an MSSM framework. However, 
we also note that these corrections introduce a huge new set of parameters into charged Higgs searches, where it is not even clear 
if charged Higgs searches at the LHC should be part of the search for supersymmetry or something more general. We leave it to 
the reader to decide this issue, and present the size of the additional supersymmetric corrections in this section. 

Just like in mass renormalization, the supersymmetric external bubbles can be resummed. It has been shown that this resum¬ 
mation consistently takes into account the leading contributions in tan/3 ll^ . Note that in the large tan /3 limit we approximate 
the off-diagonal matrix element as —rubpian(3 and it is an unsolved problem how to treat the iribAb term consistently. We 
will not resum this contribution. The difference is in any case numerically small, because we are interested only in the region 
of parameter space where tan/3 is large enough to also ensure /rtan/3 ^ Ab- If we start from the fixed relation between the 
bottom quark mass and the bottom Yukawa coupling (which we of course break in our calculation with a vanishing bottom mass 
and a running MS bottom Yukawa coupling) we can regard the resummed bottom self-energy diagrams as a shift in the Yukawa 
coupling as compared to the mass or its hxed order value: 
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The function B{p^, mi, m 2 ) is the usual scalar two-point function; Cf = 4/3 is the Casimir factor in the fundamental repre¬ 
sentation of SU{3). As mentioned above, there are similar additional terms proportional to the strong coupling or to the top 
quark Yukawa coupling, but the Af, correction is the leading contribution for large tan /3. We should also note that the resummed 
definition of these corrections is well defined for all Af, > —1, whereas the strictly fixed order series which yields (1 — 4Af,) 
breaks down for Af, > 1/4. 

We show the numerical results for the supersymmetric corrections to the total rate in Tab. |I] To illustrate the effect of the 
supersymmetric corrections we follow the SPS parameter points, which are obviously not designed to advertise searches for 
charged Higgs pairs in bottom-parton scattering. Apart from the fact that almost all 2HDM production cross sections are very 
small, we indeed see that the Af, corrections are by far dominant, and that the remaining explicit diagrams contribute at a level 
well below the total cross section theoretical uncertainty. The fact that all Af, corrections are negative reflects the fact that all 
SPS points are chosen with p > 0, a preference triggered by indirect constraints, such as measurements of the muon anomalous 
magnetic moment [g^ — 2) or the rare decay b —> 57 . 

Of course, all of the processes evaluated in section lTVI receive explicit supersymmetric corrections. For the Drell-Yan process 
we know that their effects are small Ooll . because there are no couplings enhanced by tan/3. For gluon fusion, the effects of 
squark loops can be very large if the squarks are light. We do not include these additional diagrams in this paper, because they 
require a careful analysis of the MSSM parameter space. Instead, we refer the reader to dedicated studies in Refs. EQ. We 
should also mention that fairly light charged Higgs bosons in an MSSM framework can be visible in cascade decays of heavy 
squarks and gluinos Because the production rates of squarks and gluinos are of the order of 0(10 pb) at the LHC, even a 
small branching fraction to charged Higgs bosons can be very promising. 


VI. CONCLUSIONS 

To obtain a reliable prediction for the rate of charged Higgs boson pair production at the LHC we have computed the process 
bb —> H^H~ to NLO. In previous studies E, the bottom-parton scatteringprocess had been the source of large theoretical 
uncertainty. With our improved understanding of bottom parton densities @0 1^ . we now have a recipe for calculating 
processes involving initial-state bottom quarks. To test the perturbative behavior of the LO prediction for bb H^H~, we 
study NLO corrections to the total rate as well as to the kinematical distributions of the heavy hnal state Higgs bosons. Moreover, 
we use the NLO calculation to carefully verify the validity of the bottom-parton approach. 

With this reliable prediction for the total cross section of the bottom-induced production process at hand, we compare its rate 
to Drell-Yan pair production via an s-channel gauge boson. This process is available at NLO through Prospino2 The third 
production mechanism which we include in the comparison is the gluon-fusion process with a one-loop amplitude including 
a mixed top-bottom loop E- Like the bottom-parton scattering process its rate grows like tan/3‘^ for large values of tan/3. 
We find that for low values of tan /3 the Drell-Yan process is by far dominant, and that for large values of tan /3 the loop- 
induced process overwhelms the other two. Bottom-induced production will gain relative importance only if supersymmetric 
loop contributions reduce the rate of the top-bottom loop II^IT^ . or if the bottom-parton induced rate is enhanced by universal 
supersymmetric corrections to the bottom Yukawa coupling. Finally, we show how a combined sample of charged Higgs pair 
events would look at the LHC, including the distinctly different transverse momentum spectra of the three different production 
processes. 

Note added: While we were in the process of hnalizing this work, another paper 1^ appeared in which the authors present 
results which correspond to parts of our sectionslIIlandIVl 
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